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We present a study of the proximity effect and the inverse proximity effect in a superconductor|ferromagnet 
bilayer, taking into account several important factors which mostly have been ignored in the literature so far. 
These include spin-dependent interfacial phase shifts (spin-DIPS) and inhomogeneous textures of the magneti- 
zation in the ferromagnetic layer, both of which are expected to be present in real experimental samples. Our 
approach is numerical, allowing us to access the full proximity effect regime. In Part I of this work, we study 
the superconducting proximity effect and the resulting local density of states in an inhomogeneous ferromagnet 
with a non-trivial magnetic texture. Our two main results in Part I are a study of how Bloch and Neel domain 
walls affect the proximity-induced superconducting correlations and a study of the superconducting proximity 
effect in a conical ferromagnet. The latter topic should be relevant for the ferromagnet Ho, which was recently 
used in an experiment to demonstrate the possibility to generate and sustain long-range triplet superconducting 
correlations. In Part II of this work, we investigate the inverse proximity effect with emphasis on the induced 
magnetization in the superconducting region as a result of the "leakage" from the ferromagnetic region. It is 
shown that the presence of spin-DIPS modify conclusions obtained previously in the literature with regard to the 
induced magnetization in the superconducting region. In particular, we find that the spin-DIPS can trigger an 
anti-screening effect of the magnetization, leading to an induced magnetization in the superconducting region 
with the same sign as in the proximity ferromagnet. 

PACS numbers: 74.20.Rp, 74.50,+r, 74.20.-z 



I. INTRODUCTION 

The interplay between ferromagnetism and superconduc- 
tivity has over the past decade attracted much interest from 
the condensed-matter physics community. Research on 
superconductor | ferromagnet (S|F) heterostructures continues 
to benefit from great interest, which is fueled by the exciting 
phenomena arising from a fundamental physics point of view 
in addition to the prospect of harvesting functional devices in 
low-temperature nanotechnology. 

There is currently intense activity in this particular re- 
search area (see e.g. Refs. 1,2 and references therein). 
The interest in S|F hybrid structures was boosted at the 
beginning of this millenium, primarily due to the the- 
oretical proposition of proximity-induced odd-frequency 
correlations 3 and the experimental observation of 0-ir os- 
cillations in S|F|S Josephson junctions. 4 A large amount 
of work has been devoted to odd-frequency pairing (see 

^^5.6.7.8,9,10,11,12,13,14.15,16,17.18,19,20,21) an( j ^ physics of Q _ n 
Oscillations (see e ,g, 22.23,24,25.26.27,28,29,30,31,32.33,34,35,36,37,38) j n 

S|F heterostructures. The concept of odd-frequency pairing 
dates back to Refs i 39 ' 40 ' 41 ' 42 and was recently re-examined in 
RefA 

So far, the proximity effect has received much more atten- 
tion than the inverse proximity effect. In S|F bilayers, the 
proximity effect causes superconducting correlations to pen- 
etrate into the ferromagnetic region 1 . Similarly, the inverse 
proximity effect induces ferromagnetic correlations in the su- 
perconducting region near the interface region /' 45 ' 46 ' 47 Of- 
ten, the bulk solution is employed in the superconducting re- 
gion, such that both the induced magnetic correlations and the 
self-consistency of the superconducting order parameter are 



neglected. However, it was shown in Ref. that the induc- 
tion of an odd-frequency triplet component would lead to a 
finite magnetization in the superconducting region close to 
the S|F interface. Prior to this finding, some experimental 
groups had reported findings which pointed towards precisely 
such a phenomenon 64 ' 65 . Very recently, Xia et a/. 63 presented 
an experimental observation of the inverse proximity effect 
in Al/(Co-Pd) and Pd/Ni bilayers by measuring the magneto- 
optical Kerr effect. Their data could be roughly fitted to the 
predictions of Ref. 62 , and other experiment s 44,45 ' 64 ' 65 ' 66 have 
also addressed aspects of the inverse proximity effect S|F bi- 
layers. 

In Ref. 67 , the authors investigated the proximity-induced 
magnetization in the superconducting region of a S|F bilayer, 
and found that the magnetization would oscillate in the clean 
limit (see also Ref. 68 ) and decay monotonously in the diffu- 
sive limit, with a sign opposite to the magnetization in the 
bulk of the ferromagnet. The reason for this screening be- 
havior in the superconductor was attributed to a scenario in 
which the spin-f electron of a Cooper pair near the interface 
would prefer to be located in the ferromagnetic region, while 
its spin- 1 partner would remain in the superconducting region, 
thus creating a magnetization with an opposite sign compared 
to the ferromagnet. By considering the weak proximity effect 
regime in the diffusive limit, both Ref. 62 and Ref. 67 arrived 
at this conclusion. However, it would be desirable to go be- 
yond the approximation of a weak proximity effect employed 
in previous work to investigate if this may alter how the in- 
duced magnetization in the superconducting region behaves. 

Moreover, none of the above works on the inverse proxim- 
ity effect have properly included an important property which 
is intrinsic to S|F interfaces, namely the spin-dependent in- 
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terfacial phase shifts (spin-DIPS) that occur at the interface. 
The spin-DIPS have been shown to exert an important influ- 
ence on various experimentally observable quantities in S|F 
bilayer s 28 ' 48 - 69 , and should be taken into account. For instance, 
the anomalous double peak structure in the local density of 
states (LDOS) in a diffusive S |F bilayer reported very recently 
by SanGiorgio et aim Ref. 50 was reproduced theoretically in 
Refi^ 9 - by using a numerical solution of the Usadel equation 
when including the effect of the spin-DIPS. 

So far, due to the complexity of the problem, several as- 
sumptions have been usually made when treating S|F hybrid 
structures. For instance, since the quasiclassical equations 
become quite complicated for inhomogeneous ferromagnets, 
they have been linearized in most of the previous works. How- 
ever, presently, the direction of this research field tends to- 
wards a more realistic description of S|F structures than the 
simplified models that mostly have been employed up to now. 
It is obvious that this is a necessary step in order to reconcile 
the theoretical predictions with experimentally observed data. 

Our motivation for this work is to examine the effect of in- 
homogeneous magnetization textures and spin-DIPS on both 
the proximity effect and the inverse proximity effect in S|F 
bilayers. This is directly relevant to two recent experimental 
studies 57,63 which studied the superconducting proximity ef- 
fect in the conical ferromagnet Ho and the inverse proximity 
effect in the superconducting region of a S|F bilayer, respec- 
tively. As we shall show in this work, non-trivial magnetiza- 
tion textures and spin-DIPS have profound influence on the 
physical properties of S|F bilayers, suggesting that their role 
must be taken seriously. 

We divide this work into two parts which are devoted to the 
proximity effect in the ferromagnetic region (Part I) and the 
inverse proximity effect in the superconducting region (Part 
II). In Part I, we present results where we treat the role of mag- 
netic properties at the interface and the possibility of inhomo- 
geneous magnetization thoroughly. We study the proximity- 
induced density of states (DOS) in a S|F bilayer which takes 
into account the presence of spin-DIPS at the interface and 
also the possibility of having a non-trivial magnetization tex- 
ture (such as a domain wall) in the ferromagnetic region. In 
order to access the full proximity effect regime, we do not re- 
strict ourselves to any limiting cases. Rather, we employ a 
full numerical solution of the DOS by means of the quasiclas- 
sical theory of superconductivity. We apply our theory to two 
cases of ferromagnets with an inhomogeneous magnetic tex- 
ture, namely on one hand ferromagnets with domain walls and 
on the other hand conical ferromagnets. 

In Part II, we study numerically and self-consistently the 
inverse proximity effect in a S|F bilayer of finite size upon 
taking properly into account the spin-DIPS that occur at the 
S|F interface. Our main objective is to study the influence 
exerted on the inverse proximity effect by the spin-DIPS. Sur- 
prisingly, we find that the spin-DIPS may invert the sign of 
the proximity-induced magnetization in the superconducting 
layer compared to the predictions of Refs . 62 i 67 . Consequently, 
the spin-DIPS can trigger an anti-screening effect of the mag- 
netization, which suggests that their role must be taken seri- 
ously in any attempt to construct a theory for the inverse prox- 



imity effect in S|F bilayers. We also explain the basic mecha- 
nism behind the sign-inversion induced by the spin-DIPS. 

This paper is organized as follows. In Section III Al we 
present the theoretical framework we use to perform our com- 
putations in Part I, namely the quasiclassical theory of super- 
conductivity in the diffusive limit for an inhomogeneous fer- 
romagnet using the Ricatti parametrization. In Section III Bl 
we present our numerical results for proximity-effect and the 
local density of states for the two cases of ferromagnets with 
domain walls and with conical magnetic textures. In Section 
III CI we present a discussion of our results obtained in Part I. 
Moving on to Part II of this work, we introduce a slightly dif- 
ferent notation and parametrization for the Green's function 
in Sec. IIII Al which is easier to implement for a homoge- 
neous S|F bilayer. In Sec. IIII B I we present our results for 
the inverse proximity effect, manifested through an induced 
magnetization in the superconducting region and in particular 
how it is influenced by the presence of spin-DIPS. The results 
for Part II are discussed in Sec. IIII CI and we conclude with 
final remarks in Sec. [IV] Throughout the paper, we will use 
boldface notation for 3-vectors, . ? . for 4x4 matrices, and 
for 2x2 matrices. 



II. PROXIMITY EFFECT IN A S | F BILAYER WITH AN 
INHOMOGENEOUS MAGNETIZATION TEXTURE 

A. Theoretical framework 

In the first part of our work, we shall consider the proxim- 
ity effect in the ferromagnetic region of an S|F bilayer when 
the magnetization texture is inhomogeneous. This is the case 
e.g. in the presence of a domain-wall structure or conical fer- 
romagnetism, which both will be treated below. We will use 
the quasiclassical theory of superconductivity 51 , and consider 
the diffusive limit described by the Usadel equation 5 ^. 

1. Quasiclassical theory and Green 's functions 

To account for an inhomogeneous magnetization in the fer- 
romagnet, it is convenient to parametrize the Green's function 
to obtain a simpler set of equations to solve. One possibility 
is to use a generalized ^-parametrization 53 , as follows 

A _ / M cao + (M ■ a)s p+ \ 

9 ~\ p- -M Q cao-(M-a)*sJ ' 

= c[l(M z 02 - Myas) ± M x oq] ± M £lS, (1) 

where cr,- are the identity (j = 0) and Pauli (j = 1,2,3) 
matrices, and 

3L = (£1,02, £3)- (2) 

Also, s = sinh(6>) and c = cosh((9). The Green's function is 
then completely determined by the complex functions 9, AIq, 
and M with the additional constraint M§ — M 2 = 1 in order 
to satisfy g 2 = 1, 
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However, for our purpose we find it both more convenient 
and elegant to use a Ricatti-parametrization of the Green's 
function as follows 54,55 , as follows 

/&C1-77) 2=^7 \ m 

This parametrization facilitates the numerical computations, 
and also ensures that g 2 = 1. The unknown functions 7 and 
77 are key elements in this parametrization of the Green's 
function, and will be solved for below. Here, denotes a 
2x2 matrix and 

M= (l+Tir 1 ^ (1+77T 1 - (4) 

In order to calculate the Green's function g, we need to 
solve the Usadel equation with appropriate boundary condi- 
tions at x = and x = (If- The two natural length scales 
associated with each of the long-range orders are the super- 
conducting and ferromagnetic coherence lengths 

£5 - tJDs/Ao, £f = VDp/ho, (5) 

where Ao and ho denote the bulk values of the gap and the 
exchange field. We set Dp — Dg — D for simplicity. The 
Usadel equation reads 

Dd(gdg) + i[ep 3 + diag[h ■ <r, (h ■ <r) r ] ,g]=0, (6) 

and is supplemented with the boundary condition s 28 ' 48 

2(d F gdg = [<?bcs,#] + i(G0/G T )[diag(T3,T3_),g] (7) 

at x = where the interface is spin polarized along the z-axis, 
and gdg = at x = dp. Here, d = and we define 

C = Rb/Rf (8) 

as the ratio between the resistance of the barrier region and the 
resistance in the ferromagnetic film (note that Rb = G^ 1 ). 
The barrier conductance is given by^ 8 - 

JV 

Gt^GqY^Tu, (9) 

n 

where Gq — e/h and T n is the transmission coefficient for 
channel n. The boundary conditions Eqs. ( l22l and d23l are 
derived under the assumption that T n <C 1, but this does not 
necessarily mean that the barrier conductance is small since 
there may be a large total number of channels N through 
which transport may take place. The parameter G^ describes 
the spin-DIPS taking place at the F side of the interface. 56 
Since its exact value depend on the microscopic properties of 
the barrier region, they are here treated phenomenologically. 
We finally underline that the boundary conditions above are 
valid for planar diffusive contacts. 

Since we employ a numerical solution, we have access to 
study the full proximity effect regime and also an, in principle, 
arbitrary spatial modulation h = h(x) of the exchange field. 



This is desirable in order to clarify effects associated with non- 
uniform ferromagnets, such as spiral magnetic ordering or the 
presence of domain walls. Inserting Eq. (O into Eq. (Bit , we 
obtain the transport equation for the unknown function 7 (and 
hence 77) 

D[d 2 2 + {dy)£{ d i)] + i[2e7 + h ■ (gry — ja*)] = 0, 

(10) 

with T_ = — 27V-7. The boundary condition at x = reads 

2(d F d x 2 = [2c7 - siT2 + 2( sl l2.)l} 

+ i(G /Gt)(t37-7I3), (ID 

while d x j = at x = d. For 7, we obtain 

D[d 2 j + (df)£.(df)} + i[2£7 + h ■ (7<x - g^f)} = 0, 

(12) 

with the corresponding boundary condition 

2(d F d x 2 = [2c7 - siT2 + 7(5122)7] 

-i(G /G T )(737-7 Z 3). (13) 

We have defined T_ = — 2Afi. Note that we use the bulk 
solution in the superconducting region, which is a good ap- 
proximation when assuming that the superconducting region 
is much less disordered than the ferromagnet and when the in- 
terface transparency is small, as considered here (see detailed 
discussion in Sec. IIIQ . One finds that 

-y — Zj — ( s/(l + c)\ 

The normalized DOS is finally evaluated by 

N(e)/N = Tr{Re^(l - 77)]}/2- (15) 

In what follows, we will omit the effect of spin-flip and spin- 
orbit scattering to reduce the number of parameters in the 
problem. In comparison with real experimental data, how- 
ever, the effects of these pair-breaking mechanisms are easily 
included in our framework by adding two terms cr s f and a so 
in Eq. (f2Tb (see e.g. Ref. 20 for a detailed treatment). In this 
paper, we will focus on the role of the phase-shifts obtained at 
the interface due to the spin-split bands and the inhomogene- 
ity of the exchange field in the ferromagnet. 

2. Inhomogeneous magnetization 

We will consider three types of inhomogeneous magnetic 
structures: Bloch walls, Neel walls, and conical ferromagnets 
(see Fig. Q]). An example of the latter is the rare-earth heavy 
fermion elemental magnet Ho, although we hasten to add that 
while Ho features strong ferromagnet, we will consider the 
weakly ferromanetic case. These structures are shown in Fig. 
[T]and are to be contrasted with the usual assumption of a ho- 
mogeneous exchange field in the ferromagnetic region. For 
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FIG. 1: (Color online) The three types of inhomogeneous ferromag- 
nets we will consider in this work: Bloch walls, Neel walls, and 
conical ferromagnets (such as Ho). 



the first two cases, the domain wall has a width dw and is 
taken to be located at the center of the ferromagnetic region 
(x = (If/2). The Bloch wall is thus modelled by 

h = h(cos9y + sinflz), (16) 

while y — > x for the Neel wall. Here, we have defined 

6 = -arctan[(a; - dp/2)/dw], (17) 

similarly to Refi^ 5 -. 

In the case of a conical ferromagnet, cf. Fig. Q] the mag- 
netic moment belongs to a cone. In Ho, the opening angle 
is a — Att/9 and the magnetic moment then rotates like a 
helix along the c-axis with a turning angle 9 — tt/6 per in- 
teratomic layer with distance a (see RefiS for a further dis- 
cussion). Above 21 K, the conical ferromagnetic structure 
transforms into a spiral antiferromagnetic structure. Instead of 
using an abrupt change in the magnetization direction at each 
interatomic layer, we will model this transition smoothly since 
the effective field felt between the layers probably should be 
a weighed superposition of the exchange fields from the two 
closest layers. In the ferromagnetic phase, the spatial variation 



h = h 



cos ax + sin a 



y + «»(^)*} 



(18) 



B. Results 

In what follows, we will choose the parameters of our 
model, corresponding to a realistic experimental setup in or- 
der to make our study directly relevant for experiments on S |F 
bilayers. The numerical treatment makes use of built-in rou- 
tines in MATLAB for a two-point boundary value problem 
for an ordinary differential equation. More specifically, we 
use a finite difference code which implements a three-stage 
Lobatto-Illa formula. An initial guess for the Ricatti-matrices 
is supplied with fixed boundary conditions, and the Usadel 
equation is then solved in the entire ferromagnetic region. 

In the first part of this section, we will study the effect of 
domain walls in weak ferromagnets. Weak ferromagnetic al- 
loys such as PdNi or CuNi are commonly employed in ex- 
periments, and the corresponding exchange field h depends 
on the concentration of Ni, reaching up to tens of meV. The 
modification of the DOS is most dramatic in the case when 
the energy scales for the superconductivity and the ferromag- 
netism are of the same order, h ~ A. This scenario appears to 
have been realized in Ref. 4 where Qii_ x Ni x with x = 0.44 
was used. The diffusion constant in the weakly ferromagnetic 
alloys is usually of order D ~ 10~ 4 m 2 /s. The superconduct- 
ing region is considered to act as a reservoir with thickness 
ds ^ £s> while we fix the thickness of the ferromagnetic 
region at dp/£,s = 0.5. This typically corresponds to a thick- 
ness of the ferromagnetic layer ~ 10 nm. The remaining pa- 
rameters are then the domain wall thickness dw and the term 
G$ accounting for the spin-dependent phase-shifts at the inter- 
face. Below, we will contrast a thin domain wall (dw <C c^f) 
with a thick domain wall (dw — dp) and investigate the role 
of C?0. In what follows, we choose ( = 5 corresponding to a 
situation where Rb 3> Rf- 

In the second part of this section, we will study coni- 
cal ferromagnetism, of a similar kind to that realized in the 
heavy rare earth element Holmium (Ho) under certain con- 
ditions. Recently, it was strongly suggested by experimen- 
tal data that a long-range triplet superconducting component 
was generated and sustained in a superconductor|Ho proxim- 
ity structure 57 . The experimental samples used in RefjS did 
not appear to fall into the diffusive motion regime, since Ho is 
a strong ferromagnet. More specifically, it was estimated that 
hr ~ 10 in Ref. 57 , suggesting that one would have to revert 
to the more general Eilenberger equation in order to study the 
proximity effect in Ho. In this work, we will study a coni- 
cal ferromagnet under the assumption that the diffusive limit 
is reached. For the actual structure of the magnetization, we 
choose the same parameters for Ho as those reported in Ref. 57 : 
a = 4tt/9, 6 = tt/6, and a = 0.526 nm (see Fig. [B- How- 
ever, we choose the exchange field much weaker than in Ho, in 
order to justify the Usadel approach. Thus, our results may not 
be directly applicable to Ho. While in Ref. 57 it was estimated 
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that h ~ 1 eV, corresponding to an exchange field compara- 
ble in magnitude with the Fermi energy, we choose h/A = 5 
in our study of conical ferromagnetism to ensure the validity 
of the quasiclassical approach. Assuming that £,5 = 20 nm, 
which should be reasonable for a moderately disordered con- 
ventional superconductor, we obtain a/£s = 0.0263. 

1. Domain wall 

Before proceeding to a dissemination of our results, it 
should be noted that we find identical results for the Bloch 
and Neel wall cases. This seems reasonable, since the only 
difference between those two cases is that the y-component of 
the magnetization is exchanged with the x-component. The 
long-range triplet component comes about as long as only one 
of these is non-zero, and it does not matter which one it is. 
It is also necessary for the magnetization to vary directionally 
with the ^-coordinate in order to generate the inhomogeneity 
required for the long-range triplet component. Note that the z- 
component of the magnetization is the same for the Bloch and 
Neel walls. In what follows, we only consider the Bloch wall 
configuration since the results for the Neel wall are identical. 
We also note that in our study, the magnetization is always in- 
homogeneous in the direction perpendicular to the interface, 
i.e. upon penetrating into the ferromagnetic region. In the 
case where the inhomogeneity of the magnetization is in the 
transverse direction (parallell to the interface), i.e. there is 
no variation in the x-direction, the proximity effect does not 
become long-ranged even if equal-spin correlations may be 
generated 35 . The general condition for a long-range proxim- 
ity effect is that there exists a misalignment between the triplet 
anomalous Green's function vector and the exchange field. 



We first study the thin-domain wall case d\y /dp = 0.2. To 
begin with, we shall consider the energy-resolved DOS in the 
center of the domain wall (x = dp/2) for several values of 
the exchange field. This is shown in Fig. [2] As seen, the 
zero-energy DOS is enhanced in all cases due to the pres- 
ence of odd-frequency correlations . 8 ' 9 ' 12 ' 58 The influence of 
the spin-DIPS (G$) seems to be an induction of additional 
peak features in the subgap regime. This effect is most pro- 
nounced at low exchange fields (in particular h/ Aq — 0.5 in 
Fig. [2). A possible physical explanation for the additional 
peak features in the LDOS may be the fact that acts as 
an effective exchange field in both the superconducting and 
ferromagnetic layers. 48 It thus conspires with the intrinsically 
existing exchange field in the ferromagnetic layer to yield a 
modified value of the total exchange field. This explanation is 
consistent with the fact that the position of the peaks change 
upon increasing G^. More specifically, the spin-DIPS appear 
to enhance the exchange field since the peaks move outwards 
toward the gap edge. 

Next, we investigate the thick domain wall case, and choose 
dw/d>F = 0.8. In Fig. [3j we again consider the energy- 
resolved LDOS in the middle of the ferromagnetic layer 
(x/dF = 0.5) for three different values of the exchange field. 
Upon comparison with Fig. [2] it is seen that the general trend 
upon increasing the domain wall thickness is an overall en- 
hancement of the proximity effect. The qualitative features in 
Fig. [3] are very similar to those in the thin domain wall case, 
but the enhancement at zero-energy tends to be larger, partic- 
ularly so for large values of h/Ao- Again, it is seen that the 
effect of the spin-DIPS is a modification of the total exchange 
field, amounting to a double-peak structure at subgap energies 
in the LDOS. 




FIG. 2: (Color online) Plot of the energy-resolved LDOS evaluated at x/d,F = 0.5 in the case of a thin domain-wall dw/dp = 0.2. We 
consider three values of the exchange field h and also investigate how the LDOS changes with the phase-shift G$ at the interface. 



It is also interesting to consider the spatial dependence of the zero-energy DOS in the ferromagnetic region. By using 
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FIG. 3: (Color online) Plot of the energy-resolved LDOS evaluated at x/d,F = 0.5 in the case of a thick domain-wall, dw /<1f = 0.8. We 
consider three values of the exchange field h and also investigate how the LDOS changes with the phase-shift G^ at the interface. 



local STM-techniques, it is possible to probe the DOS at (in 
principle) any location in the ferromagnetic film. The spe- 
cific choice of e = is particularly interesting in terms of 
the DOS, since it is strongly influenced by the presence of 
odd-frequency correlations. As pointed out in Refs. [T2II20I 
the behaviour of the DOS at e = may be interpreted as a 
competition between spin-singlet even-frequency correlation 
and spin-triplet odd-frequency correlation. The former tend 
to give a minigap in the DOS for subgap energies, while the 
latter yields a zero-energy peak in the DOS. Clearly, these two 
effects are competing with each other since they have a de- 
structive interplay. In the present case, one would expect that 
the domain wall structure should favor the generation of the 
odd-frequency triplet components, thus enhancing the LDOS. 
This conjecture is supported by Figs. |2]and[3] 

In Fig. 2J we plot the spatially-resolved LDOS at e = for 
several values of dw to probe directly how the odd-frequency 
correlations are affected by the domain wall thickness. As 
compared to Figs. [2] and [3] we normalize the LDOS on its 
value at x = in Fig. [4] for easier comparison between dif- 
ferent values of dw, and choose = 0. From the plot, it 
is clear that the thicker the domain wall, the more strongly 
enhanced the zero-energy DOS. This also supports the notion 
that the magnetically inhomogeneous structure favors the gen- 
eration of odd-frequency triplet components. The concomi- 
tant enhancement of the DOS may then be seen at increasingly 
larger penetration depths in the ferromagnet when the domain 
wall thickness is increased. 



2. Conical ferromagnetism 

We now turn to a study of how the superconducting proxim- 
ity effect is manifested in a ferromagnet with a conical mag- 
netization such as Ho. We fix the exchange field at h/Ao = 5 
and study how the DOS changes upon increasing the ferro- 
magnetic layer thickness. The motivation for this is to obtain 
a better understanding of how the DOS changes when only 
the long-range triplet components are present in the sample. 
In an inhomogeneous ferromagnet, the singlet component and 
the S z = triplet component are short-ranged, and pene- 
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FIG. 4: (Color online) Plot of the zero-energy LDOS induced in the 
ferromagnet for h/Ao = 5 and dp/£,s = 0.5. The lines correspond 
to dw/dF in the range [0.1,0.9] in steps of 0.1 along the arrow. 
Here, G^ is set to zero. 

trate in a distance £p = y/D/h into the ferromagnet. The 
S z = ±1 triplet components, however, are not subject to the 
pair-breaking effect originating with the Zeeman splitting, and 
can thus penetrate a much longer distance £jv = y/D/T into 
the ferromagnet, where T is temperature. Therefore, by mak- 
ing the ferromagnetic layer thick enough, one can be certain 
that there is no contribution from either the singlet or S z = 
triplet components. Since we have chosen h/Ao = 5, we find 
that the penetration depth of these components in the ferro- 
magnetic layer should be 0.44£s. 

We next turn to a study of the proximity-induced LDOS. 
In Fig. [5] we plot the energy-resolved LDOS for three layer 
thicknesses: i) d/£s = 0.1, ii) d/Cs = 0.5, and Hi) d/£s — 
0.9. In case i), both short-ranged and long-ranged compo- 
nents should contribute significantly to the LDOS. In case 
ii), the long-ranged components should dominate over the 
short-ranged ones, while finally in case Hi) only long-ranged 
components remain. This is because we evaluate the energy- 
resolved DOS at the F|I interface, x — dp, as was also done 
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in the experiment of Refs . " . 

As seen in case ii) and Hi), a pronounced zero energy peak 
is present, bearing witness of the odd-frequency correlations 
in the system. The peak is more pronounced with increasing 
thickness, since the long-range triplet correlations dominate 
over the even-frequency singlet Green's function as the thick- 
ness increases. However, case i) is qualitatively different from 



the two other thicknesses. In this case, the low-energy LDOS 
is completely suppressed in the regime G^/Gt < 1, and sud- 
denly reappears for G$ / Gt > 1. It is very interesting to note 
that the same effect was recently discovered for an S|N junc- 
tion with a magnetically active interface 60 , but in that case the 
effect was completely independent of the junction thickness. 




FIG. 5: (Color online) Plot of the LDOS at x = cLf for a conical ferromagnet with h/Ao = 5 for several values of the ferromagnetic layer 
thickness cLf- In each case, we investigate the role of the spin-DIPS (G^) at the S|F interface. 



In order to investigate this effect further, we focus on the 
zero-energy LDOS in the thin junction case in Fig. [6] As seen, 
for sufficiently thin layers dp/^s «C 1, an abrupt crossover 
takes place at a critical value of G^/Gt, qualitatively altering 
the LDOS at zero-energy. Remarkably, we find that a sim- 
ilar transition takes place upon increasing the ferromagnetic 
layer thickness. Consider a plot of the zero-energy LDOS 
in Fig. |7] as a function of dp/£s- As seen, at a critical 
layer thickness, the zero-energy LDOS rises abruptly from 
zero and acquires the usual oscillating behavior. To see how 
the full energy-resolved LDOS evolves with increasing G^ 
for a fixed thickness dp/£s = 0.1, consider Fig. [8] As 
seen, the LDOS changes qualitatively above a critical value 
of G^/Gt ~ 1.14. 

To summarize the findings of Figs. [6] |7] and [8] we have 
found that there is an abrupt crossover from a fully suppressed 
LDOS to a finite LDOS which appears at a critical thickness 
of the ferromagnetic layer, and the particular value of the crit- 
ical thickness depends on the value of G^. In a similar way, 
we find that there is an abrupt change appearing at a critical 
value of Gcf, for sufficiently thin layers. The natural question 
is: what is the reason for these changes? An important clue is 
found in the fact that when the LDOS is fully suppressed, the 
odd-frequency correlations must be zero 12 . The presence of 
odd-frequency correlations will in general lead to an enhance- 
ment of the LDOS at zero-energy, which at present is one of 
the main suggestions put forth in the literature with regard to 
the issue of how to obtain clear experimental signatures of 
this exotic type of superconducting pairing. Therefore, the 
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FIG. 6: (Color online) Plot of the zero-energy LDOS at x = dF for 
a conical ferromagnet with h/A — 5 as a function of the normal- 
ized spin-DIPS parameter G^/Gt- Below a critical value for dF, a 
qualitatively different behavior is observed for the LDOS. 



abrupt transition from a fully suppressed LDOS to a LDOS 
which is enhanced even compared to the normal-state value 
is a strong indicator of a symmetry-transition from the usual 
even-frequency correlations to a state of mixed even- and odd- 
frequency correlations, or possibly even pure odd-frequency 
correlations. It is therefore clear that the spin-DIPS occuring 
at the interface have paramount consequences with regard to 
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FIG. 7: (Color online) Plot of the zero-energy LDOS at x = (If as 
a function of d^/Cs for several values of G$. As seen, there is a 
critical thickness at which the zero-energy LDOS becomes non-zero. 
We have used h/Ao — 5. 



the symmetry-properties of the induced superconducting cor- 
relations in the ferromagnet. Due to the complexity of the 
problem, it is unfortunately not possible to give an exact an- 
alytical treatment of the influence of on the symmetry- 
properties of the anomalous Green's function. 
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depends on the structure of the magnetic texture, which is de- 
termined by the parameters {a, a, 9} in Fig. Q] We here focus 
on the role of a and 9, which control respectively the direc- 
tion and the speed of rotation of the magnetization upon en- 
tering the ferromagnetic layer. Thus, we keep a/£s fixed at 
a/£s = 0.0263. In Fig. [9] we present results for the zero- 
energy LDOS at x = dp as a function of 9 for several values 
of a. The LDOS displays oscillations as a function of 9, and 
eventually seems to sature upon increasing 9. This may be un- 
derstood microscopically by realizing that when the rotation 
of the magnetization texture becomes faster, i.e. increasing 9, 
the effective magnetization felt by the Cooper pair averages 
out to zero for the rotating components. For our setup, this 
would mean that only the h x -component should remain non- 
zero, while h y = h z = 0. To verify this scenario, we have 
also plotted the results in the h y = h z = case in Fig. [^(dot- 
ted lines) for each value of a, which is seen to coincide with 
the limiting behavior in the the high-^ case. It is interesting to 
note that for a = tt/2, the LDOS vanishes completely above 
a critical value for 9. This may be understood by noting that 
h x = when a — tt/2. Thus, when 9 increases, we have 
ih x ) = (h y ) = 0, causing the ferromagnetic layer to act as a 
normal metal. 
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FIG. 9: (Color online) Plot of the zero-energy LDOS at x — (If 
as a function of 9 for several values of a. We have fixed ft/Ao = 
5,g!f/Cs = 0.5, = to focus on the effect of the magnetic 
structure. The dotted lines give the result for h y — h z — 0, which 
corresponds to the saturating behavior when 6 increases since the 
average value of h y and h z vanishes in this limit. The inset shows 
the case a = tt/2, corresponding to h x — 0. For increasing 8, the 
ferromagnetic layer effectively acts as a normal metal, thus causing 
a complete suppression of the zero-energy LDOS. 



FIG. 8: (Color online) Plot of the LDOS at a; = d,F for a conical 
ferromagnet with h/Ao = 5 and cJf/^s = 0.1 for several values of 
the spin-DIPS (G^) at the S|F interface. In a), G^/Gt is below the 
critical value, while in b) G$ /Gt is larger than the critical value. 

In the remaining part of the discussion of conical ferromag- 
nets, we wish to focus on how the proximity-induced LDOS 



C. Discussion 

The main approximation that we have made in our calcu- 
lations is to use the bulk solution for the order parameter of 
the superconductor. Although this approximation is expected 
to be satisfactory in the regime ds S> {^Sj^f}, such that 
the superconductor acts as a reservoir, there are two aspects 
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which are lost upon doing so. One aspect is the depletion of 
the superconducting order parameter near the interface. The 
depletion may be disregarded in the tunneling limit 61 (low bar- 
rier transparency), and we do not expect that an inclusion of 
the spatial profile of the superconducting order parameter near 
the interface should have any qualitative influence upon our 
results, as long as the superconducting order parameter is not 
dramatically reduced at the interface. 

The assumption of a step-function superconducting order 
parameter is commonly employed in the literature, but let us 
for the sake of clarity here examine a bit more carefully un- 
der which circumstances this is truly warranted. In the present 
work, we have considered a superconducting reservoir of size 
ds ^ £,s an d a ferromagnetic film of size dp < £g. For a 
weak ferromagnet considered here, the ferromagnetic coher- 
ence length £f is comparable in size to £s. Also, we have 
considered the case where £ = Rb/Rf ^> 1, corresponding 
to a low barrier transparency, which should be experimentally 
relevant. To investigate quantitatively how much the super- 
conducting order parameter is suppressed near the interface, 
let us fix h/A = 10, ds/Zs = 5, d F /£, F = 1, and ( = 5. 
Using a numerical approach for S|F bilayer with a homoge- 
neous exchange field as employed in Part II of our paper, we 
obtain the gap self-consistently with the result shown in Fig. 
[TTI It is also necessary to introduce the barrier asymmetry fac- 
tor 7 = £,s&f I (S.fo's), where ctf(S) is the conductivity in the 
F (S) layer. Here, we set 7 = 1. As seen, the depletion of the 
gap is quite insensitive to the value of G^, and we have veri- 
fied that the depletion of the gap is virtually the same even up 
to ferromagnetic layer thicknesses of dp/£p — 4. As recently 
pointed out in Ref. 49 , the step-function approximation breaks 
down for low values of £ and/or high values of 7, and if the 
spin-DIPS G^ induced on the superconducting side are large 
in magnitude compared to the tunneling conductance Gt, the 
suppression of the gap becomes more pronounced. 

The second aspect which is lost is the inverse proximity ef- 
fect in the superconductor. The inverse proximity effect is, in 




FIG. 10: (Color online) Plot of the zero-energy LDOS at x = <1f in 
the case of a Bloch domain wall with h/Ao = 5. In a) and c), we 
plot the LDOS as a function of the spin-DIPS G4,, while in b) and 
d) we plot it as a function of the ferromagnetic layer thickness d F . 
For thin layers g(f/£s <C 1, one observes an abrupt transition from a 
fully suppressed DOS to a non-zero DOS at a critical value for either 
Gtj, or dis- 



similarity to the depletion of the order parameter, expected to 
be small when the interface transparency is low and ds S> dp. 
Nevertheless, the presence of the spin-DIPS at the interface, 
modelled through the parameter G^, could have some non- 
trivial impact on the correlations in the superconductor. Cot- 
tet showed that this may indeed be so in Ref.|69l at least when 
the superconducting layer is quite thin. The full effect exerted 
on the LDOS by the presence of spin-DIPS on both sides of 
the interfaces was recently investigated numerically in an S|F 
bilayer 49 . However, no study so far have investigated how the 
proximity-induced magnetization in the superconducting re- 
gion is affected by spin-DIPS. We will proceed to investigate 
this particular issue in detail in Part II of this work. 

Above, we have considered the diffusive limit £s Amp ^S> 1, 
where l lmp = vpr is the mean free path. Although the mag- 
netic texture we have considered in the second part is identical 
that of the conical ferromagnet Ho, one important difference 
is that Ho is a strong ferromagnet, contrary to the case studied 
here. This means that the diffusive limit condition hr -C 1 
is not fulfilled for Ho, and it was in fact estimated inS that 
hr ~ 10. This calls for a treatment with the more general 
Eilenberger equation, which allows for a study where the en- 
ergy scale of the Zeeman-splitting is comparable or larger than 
the self-energy associated with impurity scattering. A natu- 
ral continuation of this work would therefore be to study a 
proximity-structure of a superconductorjconical ferromagnet 
for an arbitrary ratio of the parameter hr. Such an endeavor 
would nevertheless be quite challenging unless a weak prox- 
imity effect is assumed. In the present work, we have not 
restricted ourselves to any limits with regard to the barrier 
transparency or the proximity effect. Although the exchange 
field considered for the conical ferromagnet in this paper is 
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FIG. 1 1 : (Color online) Self-consistent solution for the spatial pro- 
file of the superconducting order parameter, using the approach de- 
scribed in Part II of this paper. The choice for parameter values are 
specified in the main text. To add stability to the numerical calcu- 
lations, we added a small imaginary number \r\ to the quasiparticle 
energies e, with r\ — 0.05Ao, effectively modeling inelastic scatter- 
ing. 
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smaller than the one realized in Ho, we expect that our re- 
sults may be qualitatively relevant for STM-measurements in 
superconducting junctions with Ho. In general, increasing the 
exchange field amounts to a quantitative reduction of the mag- 
nitude of the proximity effect. 

Finally, we show that the zero-energy DOS for the do- 
main wall case exhibits a similar crossover behavior as the 
conical ferromagnetic case upon varying G^ and dp when 
^f/Cs ^ 1- I n Fig- HQ] the zero-energy DOS is plotted for 
the thick-domain wall case to illustrate this effect - the results 
are very similar even for dw/d F <C 1 when dp/^s <C 1- 
Once again, it should be noted that a complete suppression 
of the DOS amounts to pure even-frequency superconduct- 
ing correlations induced in the ferromagnetic region, since 
the presence of odd-frequency correlations enhances the zero- 
energy DOS. The exact microscopic mechanism behind the 
abrupt crossover occuring at critical values of G$ and dp, 
respectively, remains somewhat unclear. A possible resolu- 
tion to this behavior is the observation that the spin-DIPS may 
conspire with the proximity-induced minigap in the ferromag- 
netic region for sufficiently thin layers (dp/^s <C 1) and yield 
a zero-energy DOS of the form N(0) ~ l/y^ - G\, as 

noted in Ref.—. In this case, a scenario similar to the one 
of a thin-film superconductor in the presence of an in-plane 
magnetic field is realized, where the spin-resolved DOS expe- 
riences a quasiparticle energy-shift with ±h. In this case, the 
role of the exchange field is played by G^ while the role of 
the superconducting gap is played by Gp- We do not observe 
the effects shown in Fig. [I0]for larger values of dp, which is 
consistent with the fact that the minigap is completely absent 
in this regime since the proximity effect becomes weaker. 



III. INVERSE PROXIMITY EFFECT IN A S | F BIL AYER 
WITH A HOMOGENEOUS MAGNETIZATION TEXTURE 



In this part of the paper, we will consider the inverse prox- 
imity effect of an S|F bilayer, where the exchange field is 
fixed and parallel to the z-axis, manifested through an in- 
duced magnetization near the interface of the superconduct- 
ing region. We will again employ the quasiclassical theory of 
superconductivity 51 , and consider the diffusive limit described 
by the Usadel equation 52 , as this is experimentally the most 
relevant case. Our approach will be to solve the Usadel equa- 
tion and the gap equation for the superconducting order pa- 
rameter self-consistently everywhere in the system shown in 

Fig.ua 



A. Theory 

We will use the conventions and notation of Ref. 20 , which 
also allows for an inclusion of magnetic impurities and spin- 
orbit coupling if desirable. To facilitate the numerical imple- 
mentation, we employ the following parametrization of the 



Green's functions: 
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j = {S,F} (19) 



where we have introduced 

sinh(6' ( j J ), c a j — cosh^o-j). 



(20) 



Note that (gj) 2 = 1 is satisfied. The parameter 6> CTJ is a mea- 
sure of the proximity effect, and obeys the Usadel equation 

D o 9 %j + 2l ( e + sinh(6» (T , i ) 



-2iaAcosh(6» CT!i ) = 0, <t = {T,|} 



(21) 



in the superconducting (h = 0, j = S) and ferromagnetic 
(A = 0, j = F) layer, respectively. Above, D$ and Dp 
denote the diffusion constants in the two layers, e is the quasi- 
particle energy, A is the pair potential, while h is the exchange 
field. The two latter are in general subject to a depletion close 
to the S|F interface. 

The boundary condition for the ferromagnetic Green's 
function, gp, reads^ 

2£, F g F d x g F = j T [gs, 9f] + i70,f[Q!3, 9f] (22) 

at x = 0, and gpd x gp — at x = dp. Here, .?. de- 
notes a 4 x 4 matrix in spinigiparticle-hole space. Also, 
&3 = diag(l, — 1, 1, — 1). For the superconducting Green's 
function, gs, we have 

2(ts/l)gsd x gs = -Kt\9f,9s] ~ i70,s[«3,5s] (23) 

at x = 0, and gsd x gs = at x = —d$- Above, we have 
defined 

7t = G t £f/{A<t f ), 7</»,f(S) = G^^^p/iAap), (24) 

and the barrier asymmetry factor 

7 = £,svf/(£, F vs)- (25) 

Moreover, A is the tunneling contact area, while crprs) are the 
normal-state conductivities. Note that 



Aa F ( S ) 



If{S)/Rf(S) 




(26) 
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FIG. 12: (Color online) The experimental setup proposed in this pa- 
per: a superconductor | ferromagnet bilayer. 
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where dp(S) is the thickness of the layer and Rf(S) is the 
normal-state resistance. 

In total, the interface between the S and F regions is thus 
characterized by three parameters: the normalized barrier con- 
ductance 7t, the spin-DIPS 7^3 and 70,f on each side of the 
interface. In what follows, we will study the mutual influ- 
ence of superconductivity and ferromagnetism on each other, 
instead of assuming the bulk solution for gg in the supercon- 
ducting region, as is usually done in the literature. We solve 
the Usadel equation self-consistently in both the S and F layer, 
supplementing it with the gap equation: 

a ^Nf^ f d£tanh ^ £ / 2 ) yVRe{smh(0 CT )}, (27) 
2 Jo 

where we choose the weak coupling-constant and cut-off en- 
ergy to be NpX — 0.2 and uj/Aq — 75. When obtaining 
the Green's functions, a number of interesting physical quan- 
tities may be calculated. For instance, the normalized LDOS 
is obtained according to 

N(e)/N = Re{cosh6> T + cosh0;}/2. (28) 

Experimentally, the LDOS may be probed at x — —ds in the 
superconducting layer and x = dp in the ferromagnetic layer 
by performing tunneling spectroscopy through the insulating 
layer. In principle, it is also possible to obtain the LDOS at 
any position x by using spatially-resolved scanning tunneling 
microscopy. 

The quantity of interest which we shall focus on in this 
work is the proximity-induced magnetization in the supercon- 
ducting region. A few words about the sign of the magne- 
tization in the problem is appropriate. First, recall that the 
magnetic moment fi of an electron is directed opposite to its 
spin S, namely fi ~ —(e/m e )S, where e = |e and m e is the 
electron charge and mass. Therefore, if the exchange energy 
h favors spin-| electrons energetically, the resulting magneti- 
zation M of the ferromagnet will be directed in the opposite 
direction, M \\ {—z). 

In the absence of a proximity effect, we have M = in the 
superconducting region and M = MqZ in the ferromagnetic 
region, where 

M ~ -fM B N h (29) 
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FIG. 13: (Color online) (a) Plot of the proximity-induced magneti- 
zation for 7^>,F = upon varying the spin-DIPS 7^,s on the su- 
perconducting side, (b) Plot of the proximity-induced magnetization 
for 70,s = upon varying the spin-DIPS 70, f on the ferromag- 
netic side. We have used ds /£s = 0.2, and the other parameter val- 
ues are discussed and provided in the main text. Note that the lines 
with a (blue) circle marker are equal in (a) and (b), corresponding 
to 70, s = 70, f = 0. As seen, 7^,3 affects SMs/Mo much more 
strongly than j$,f- In (b), the magnetization switches sign upon 
increasing 70, f further (for the present parameters, the sign switch 
occurs around 70,f — 1-3). 



in the quasiclassical approximation h <C ep- Now, the change 
in magnetization due to the proximity effect may be calculated 
according to 

cr 

in both the superconducting and ferromagnetic region. Using 
a quasiclassical approach, the above expression translates into 
a normalized change in magnetization 

r°° de 

6M/M = - — ^crRe{cosh6» CT }tanh(/3e/2). (31) 



In the ferromagnetic region, the normalized magnetization 
M/Mq is therefore 1 + SMp/Mo, while in the supercon- 
ducting region we have an induced magnetization 5Ms/Mq, 
where SMpi s) is determined by Eq. OTb on the ferromagnetic 
(superconducting) side of the interface. 

Although we shall be concerned with a full numerical solu- 
tion when presenting our results in Sec. IIIIB1 let us for com- 
pleteness sketch how an analytical solution may be obtained 
under the assumption of a weak proximity effect. Including 
the spin-DIPS, the analytical results obtained here are thus a 
natural extension of the results in Ref. 62 , where the spin-DIPS 
were neglected. We remind the reader that spin-DIPS occur 
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whenever there is a finite spin-polarization in the ferromag- 
netic region or when the barrier itself is magnetic. 

In the weak-proximity regime, the Usadel equation in the 
ferromagnetic region becomes 



D F dl 



+ 2i(e + ah)66l = 0, 



(32) 



where the linearization of Eqs. dT9b and d2Tb amounts to 
6<j,f — * wnere \59^\ <S 1- The general solution is readily 
obtained as 



= A a {c 



(33) 



upon taking into account the vacuum boundary condition 
d x 59^ = at x = dp, and defining 

k 2 a = 2\{e + ah)/D F . (34) 

In the superconducting region, we obtain the Usadel equation 

D s dl50 s a + 2i(ec B cs - As B cs) = 0, (35) 

under the assumption that the superconducting order param- 
eter is virtually unaltered from the bulk case. This is a valid 
approximation for {7, jt} <C 1 and not too large 7^5 (typ- 
ically 70.5 < 1), which we have verified by using the full 
numerical solution. Here, 69% is the deviation from the bulk 
BCS solution, i.e. 6» CT , S -» cr9 BC s + 69% with \60%\ < 1 and 

cbcs = cosh(0 B cs), sbcs = sinh(6>Bcs), 
6*bcs = atanh(A/e). 

In this case, the general solution reads 



(36) 



-lKX — 2lK,ds 



), 



(37) 



when incorporating the vacuum boundary condition d x 89% 
at a; = —d$, upon defining 



k 2 = (ecbcs - AsbcsVAs 



(38) 



The remaining task is to determine the unknown coefficients 
{A a , B a }. Linearizing the boundary conditions Eq. d22l > and 
d23l . we obtain at x = 

£ s d x 69%/ 1 = lT {c89^ - asBCS - c BC s^) 

£ F d x 5e* = j T (cS9^ - asBes - c BCS 69%) 

+ mrfwSBZ . (39) 
From these boundary conditions, one derives that 

z% Sbcs[^3 (7tct + i70,s) - zl <*lT 



A„ = 



z 2 z 3 



os B cs7t, 



B„ 



sbcs[zi ctjt ~ z^i'jro- + 17^,5)] 



z 2 z 3 z 1 z 4 
Here, we have defined the auxiliary quantities: 



(40) 



-7TCBcs(l+e^ dp ) 



ik£s(1 - c 



— 2inds 



7 



+ c B cs(7t + 10-70,5) (1 + 1 



zl = iAv£f(1 - e 2lk ° dF ) - (7TCBCS + oi 7 0,f)(1 + 



J2\k a df. 



cbcs7t(1 + e" 



-2\K,ds 



(41) 
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FIG. 14: (Color online) Plot of the proximity-induced magnetization 
for 7^>,F = upon varying the spin-DIPS 7^s on the superconduct- 
ing side using ds/£s = 1.0. 




FIG. 15: (Color online) Plot of the spatial dependence of the total 
magnetization at zero temperature for ds/^s = 0.2 and 7<^f = 

70,s = 0.0. 



Eqs. d33l l. (l37l i. and (l40b constitute a closed analytical solution 
for the Green's functions in the entire S|F bilayer. To use this 
analytical solution, one should verify that \89^' S \ <C 1 for the 
relevant parameter regime. Spin-flip and spin-orbit scattering 
may also be accounted for in the analytical solution of the 
Green's function by adding appropriate terms to the Usadel 
equation. The calculation is then performed along the lines of 
Refs« 
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FIG. 16: (Color online) Plot of the spatial dependence of the total 
magnetization at zero temperature for ds/£s = 0.2, 70, f = 0.0 and 

70, s = 1-0. 



B. Results 

We are now in a position to evaluate the proximity-induced 
magnetization numerically. The full (non-linearized) Usadel 
equation will be employed, such that we are not restricted to 
the weak-proximity effect regime. To stabilize the numerical 
calculations, we add a small imaginary number to the quasi- 
particle energy, e — > e + vq, with rj = 0.05Ao- We focus 
on the results reported very recently by Xia et al.,— and take 
ds/£,S — 0.2, cIf/£,f — 1.0, and h/Ao — 15 as a reason- 
able set of parameters which should be relevant to this experi- 
ment. Also, we assume that the junction conductance was low, 
7t = 0.1, and we set the barrier asymmetry factor to 7 = 0.2, 
corresponding to a scenario where the superconducting region 
is much less disordered than the ferromagnetic one. We will 
also investigate the case ds/£,s — 1-0, to see how the prop- 
erties of the system changes when going away from the limit 
ds/£,S *C 1. We underline that our main objective in this 
work is to investigate the influence of the spin-DIPS on the 
proximity-induced magnetization in the system, such that we 
mainly vary 70, p and 7^,5 while keeping the other parameters 
fixed. 

Let us first consider the temperature-dependence of the 
proximity-induced magnetization in the superconducting re- 
gion in Fig. [13] To clarify the role of the spin-DIPS on 
each side of the interface, we plot 5 Ms /Mo for several val- 
ues of 70 s in Fig. Q3] (a) while keeping 70. f = fixed. 
Conversely, we plot SMs/Mq for several j^>,f m Fig- [131(b) 
with 70, s = 0. In both cases, we plot the proximity-induced 
magnetization at x = —ds- One obvious difference between 
these two scenarios is that the spin-DIPS on the superconduct- 
ing side, 70, s, influence the proximity-induced magnetization 
much stronger than 70, p. The same thing is true with regard 
to the influence of spin-DIPS on the superconducting order 
parameter: 70.5 influences the spatial profile of A much more 



than what 7^ does. From Fig. [13] it is clearly seen how 
the proximity-induced magnetization may switch sign upon 
increasing the magnitude of the spin-DIPS 70, s- We have 
checked numerically that this effect also takes upon increasing 
70,F when keeping 70, 5 = 0. Thus, increasing either 70,s or 
70. f can lead to a sign change of the proximity-induced mag- 
netization in the superconducting region. It is then clear that 
the conclusion of Ref. 67 that only spin screening is possible in 
diffusive S|F bilayers does not hold in general, since the pres- 
ence of spin-DIPS alters the screening effect. In what follows, 
we focus on the role of 70,5 since its impact on 6Ms/Mq is 
much greater than that of 70,f- In Fig- El we consider the 
case ds/£,s — 1-0 to show that the sign change of the magne- 
tization persists when going away from the limit ds/£,s *C 1. 
The spatial profile of the total magnetization in the F and S 
regions are shown for the case d$/£, — 0.2 with 70.5 = 0.0 
in Fig. [T5land 70.5 = 1.0 in Fig. [16] It is seen that the mag- 
netization decreases in a monotonic fashion toward the super- 
conducting region, and reaches its bulk value deep inside the 
ferromagnetic region. In the superconductor, magnetization is 
induced near the interface and decays with the distance from 
the interface. 



C. Discussion 

We propose the following explanation for the anti-screening 
effect observed upon increasing 70.5. The effect of the spin- 
DIPS in the case of a thin superconducting layer ds -C £s 
in Ref.® was shown to be equivalent to an internal magnetic 
exchange splitting h e g in the superconducting region. There- 
fore, the magnitude of the magnetization in the superconduc- 
tor should essentially grow with an increasing value of 70, s. 
If this is the case, the proximity-induced magnetization should 
also be sensitive to the sign of 70. 5, as the opposite spin 
species would be energetically favored when comparing the 




FIG. 17: (Color online) Plot of the proximity-induced magnetization 
at T = as a function of 7^,3 for two different values of ds/^s- 
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case 70.5 with (—70,5). To test this hypothesis, we plot in 
Fig. [17] the proximity-induced magnetization at T = as a 
function of the spin-DIPS on the superconducting side, 7^.5 
(keeping 7^ = 0). The results confirm our hypothesis - it 
is seen that SMs/Mq is an antisymmetric function of 70,s- 
The influence of j,p.s can also be seen directly in the LDOS 
in the superconducting region. For 7^.5 ^ 0, we obtain a 
double-peak structure in the LDOS at x = —ds in agreement 
with Refs i 49 ' 69 , while the superconducting order parameter de- 
pletes very little close to the interface for the chosen parameter 
values. In general, the depletion of the superconducting order 
parameter is found to be small as long as {77-, 7} < 1 and 
70 t 5 ~ 1 or smaller. 

In Ref. 62 , the inverse proximity effect of an S|F bilayer 
was studied without taking into account the presence of spin- 
DIPS, with the result that the proximity-induced magnetiza- 
tion in the superconducting region would have the opposite 
sign of the proximity ferromagnet, i.e. a screening effect. It 
was proposed in Ref. 62 that this behavior could be understood 
physically by considering the contribution to the magnetiza- 
tion from Cooper pairs which were close to the interface: the 
spin- 1 electron would prefer to be in the ferromagnetic region 
due to the exchange energy, while the spin- j electron remain- 
ing in the superconducting region then would give rise to a 
magnetization in the opposite direction of the proximity fer- 
romagnet. However, it is clear from the present study that this 
simple picture must be modified when properly considering 
the spin-DIPS 7^ 5 on the superconducting side of the junc- 
tion, since they act as an effective exchange field inside the 
superconductor. 

In this paper, we have evaluated the proximity-induced 
magnetization in the vicinity of the interface without tak- 
ing into account the Meissner response of the superconduc- 
tor. This should be permissable in a thin-film geometry as 
the one employed in Ref. 63 , where the screening currents are 
suppressed. In particular, for a field in the plane of the super- 
conducting film (see Fig. [T2| >. the Meissner effect should be 
strongly suppressed 71 for ds/£,s <C 1- 



IV. SUMMARY 

In conclusion, we have in Part I of this work investigated 
the proximity effect in a superconductor|inhomogeneous fer- 
romagnet junctions. Proper boundary conditions which take 
into account the spin dependent phase-shifts experienced by 
the reflected and transmitted quasiparticles were employed. 
As an application of our model, we have studied the LDOS in 
the ferromagnet in the presence of domain walls and a conical 
magnetic structure. We find that the presence of a domain wall 
enhances the odd-frequency correlations induced in the ferro- 
magnet, manifested through a zero-energy peak in the LDOS. 
For the conical ferromagnet, we show that the spin-dependent 
phase shifts originating with the interface have a strong qual- 
itative effect on the LDOS, especially for thin layers. In par- 
ticular, we find an abrupt crossover from a fully suppressed 
LDOS to a finite LDOS which appears at a critical thickness 
of the ferromagnetic layer, and the particular value of the criti- 



cal thickness depends on the value of G$. In a similar way, we 
find that there is an abrupt change appearing at a critical value 
of G<j, for sufficiently thin layers. We speculate that the rea- 
son for this could be a symmetry-transition from even- to odd- 
frequency correlations for the proximity-amplitudes in the fer- 
romagnetic region. The theory developed in the present paper 
takes into account both the phase-shifts acquired by scattered 
quasiparticles at the interface due to the presence of ferromag- 
netic correlations, and also an arbitrary inhomogeneity of the 
magnetic texture on the ferromagnetic side. Our results for 
the conical ferromagnetic structure should be relevant for the 
material Ho, which was used in Ref. 57 to indicate the presence 
of long-range superconducting correlations. 

In Part II of this work, we have investiged numer- 
ically and self-consistently the inverse proximity effect 
in a superconductor j ferromagnet (S|F) bilayer, manifested 
through an induced magnetization in the superconducting re- 
gion. We find that the interface properties play a crucial role 
in this context, as the spin-dependent interfacial phase-shifts 
(spin-DIPS) may invert the sign of the proximity-induced 
magnetization. This finding modifies previous conclusions 
obtained in the literature, and suggests that the influence of 
the spin-DIPS should be properly accounted for in a theory 
for the inverse proximity effect in S|F bilayers. 
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APPENDIX A: SPIN- ACTIVE BOUNDARY CONDITIONS 

To facilitate and encourage use of the spin-active boundary 
conditions required for an S|F interface, we here write down 
their explicit form in the diffusive limit for the case of a mag- 
netization in the ^-direction, following Ref i 48 i 69 . Consider a 
junction consisting of two regions 1 and 2, as shown in Fig. 
[T8l The regions have widths dj and bulk electrical resistances 
Rj . The matrices used below are 4 x 4 matrices in particle- 
holeigispin space, using a basis 

M(r,*)\ 

<rf(r,t) ■ 
Wl(r,t)/ 

Introducing a = diag(l, — 1, 1, — 1) = diag(er 3 , ct 3 ), where 
03 is the third Pauli matrix in spin-space, we may write the 
boundary conditions as follows: 

2(di/Ri)gid x gi = Gt [<7i , £2] - iG<f,,i[a, gi], 

2{d 2 J R 2 )g 2 d x g 2 = G T [gi, 92] + iG^, 2 [d, g 2 ]. (Al) 
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Here, Gt is the conductance of the junction while G<pj are 
the phase-shifts on side j of the interface. The parameters 
{Gt, G</>,j} may be calculated by relating them to micro- 
scopic transmission and reflection probabilities within e.g. a 
Blonder- Tinkham-Klapwijk (BTK)S framework. Explicitly 
spin-active barriers were considered in ballistic S|F bilayers 
using the BTK-approach for both s-wave 73 and d-wava 7 ^ su- 
perconductors. In the absence of spin-DIPS (G^j — > 0), 
Eq. ( lAU reduce to the Kupriyanov-Lukichev non-magnetic 
boundary conditions 75 . Let us make a final remark concern- 
ing the treatment of interfaces in the quasiclassical theory of 
superconductivity. We previously stated that the application 
of the present theory requires that the characteristic energies 
of various self-energies and perturbations in the system are 
much smaller than the Fermi energy £f- At first glance, this 
might seem to be irreconcilable with the presence of inter- 
faces, which represent strong perturbations varying on atomic 



length scales. However, this problem may be overcome by in- 
cluding the interfaces as boundary conditions for the Green's 
functions rather than directly as self-energies in the Usadel 
equation. 




o?l G?2 a:-axis 

FIG. 18: (Color online) lunction consisting of two regions 1 and 2 
with an interface perpendicular to the a; -axis. 
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